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\Q ' Abstract 

A methodology is presented for bounding all higher moments of the local hydrostatic stress 
field inside random two phase linear thermoelastic media undergoing macroscopic thermome- 
0^ ■ chanical loading. The method also provides a lower bound on the maximum local stress. Explicit 

formulas for the optimal lower bounds are found that are expressed in terms of the applied macro- 
. , scopic thermal and mechanical loading, coefficients of thermal expansion, elastic properties, and 

■ volume fractions. These bounds provide a means to measure load transfer across length scales 

relating the excursions of the local fields to the applied loads and the thermal stresses inside 
each phase. These bounds are shown to be the best possible in that they are attained by the 
Hashin-Shtrikman coated sphere assemblage. 

> ■ 1 Introduction 

Over the last century major strides have been made in the characterization of effective constitutive 
laws relating average fluxes to average gradients inside random heterogeneous media see for example 
O HZ1 [2H HH HI] • However much less is known about the point wise behavior of local fluxes 
and gradients fields inside random media. While it is true that efficient numerical methods capable 
of resolving local fields are available for prescribed microstructures, it is also true that for most 
applications only a partial statistical description of the microstructure is available. Thus for these 
cases one must resort to bounds or approximations for the local fields that are based upon the 
available statistical descriptors of the microgeometry and the applied macroscopic loading. Bounds 
are useful as they provide a means to quantitatively assess load transfer across length scales relating 
the excursions of the local fields to applied macroscopic loads. Moreover, they provide explicit criteria 
on the applied loads that are necessary for failure initiation inside statistically defined heterogeneous 
media pQ. In this paper we develop lower bounds on local field properties for statistically defined 
two phase microstructures when only the volume fraction of each phase is known. Here the focus 
is on lower bounds since volume constraints alone do not preclude the existence of microstructures 
with rough interfaces for which the V norms of local fields are divergent see [TH] , [3] , and also [T2] . 

We present a methodology for bounding the L p norms, 2 < p < oo, of the local hydrostatic stress 
field inside random media made up of two thermoelastic materials. The method is used to obtain 
new optimal lower bounds that are given by explicit formulas expressed in terms of the applied 
thermal and mechanical loading, coefficients of thermal expansion, elastic properties, and volume 
fractions. We show that these bounds are the best possible in that they are attained by the local 
fields inside the coated sphere assemblage originally introduced in [7]. It has been known since 1963 
that the coated spheres microstructure exhibits extreme effective elastic properties |S]. However 
it was discovered only recently in [T5], [TJ] that this geometry supports extreme local fields that 
minimize the maximum local hydrostatic field over all two phase elastic mixtures in fixed volume 
fractions. More recently several scenarios are identified for which, in the absence of thermal stresses, 
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these microstructures attain lower bounds on the total local stress field inside each material when 
the composite is subjected to mechanical loading see, pQ. 

In this paper we consider mixtures of two thermoelastic materials with shear and bulk moduli 
specified by fi±, k±, fii, fea and coefficients of thermal expansion given by hi and hi. New lower 
bounds are presented for elastically well ordered phases for which k\ > ki and \i\ > [ii as well as for 
non well ordered phases such that ki > k\ and /Ltj > /i 2 . For each of these cases we consider both 
macroscopic mechanical and thermal loads and present bounds that hold for h\ > hi and hi > h±. 
The set of bounds and optimal microstructures for the well ordered case are listed in Section 3 and 
optimal lower bounds for the non well ordered case are listed in Section 4. The methodology for 
deriving the bounds is presented in Section 5. 

The optimal bounds and the associated coated sphere microstructures given in Sections 3 and 
4 show that there are combinations of applied stress and imposed temperature change for which 
the local hydrostatic stress inside the connected phase of the coated sphere assemblage vanishes 
identically. Other loading combinations are seen to cause the stress inside the included phase of the 
coated sphere assemblage to vanish identically. Thus for these cases the applied hydrostatic stress 
is converted into a pure local shear stress inside a preselected phase. 

Recent related work provides optimal lower bounds on local fields in the absence of thermal loads. 
The work presented in [1] provides new optimal lower bounds on both the local shear stress and 
the local hydrostatic component of stress for random media subjected to a series of progressively 
more general applied macroscopic stresses. These bounds are explicit and given in terms of volume 
fractions, elastic constants of each phase, and the applied macroscopic stress. Earlier work considers 
random two phase elastic composites subject to imposed macroscopic hydrostatic stress and strain 
see, [14] and [15], as well as dielectric composites subjected applied constant electric fields see, 
[13j . Those efforts deliver optimal lower bounds on the L p norms for the hydrostatic components 
of local stress and strain fields as well as the magnitude of the local electric field for all p in the 
range 2 < p < 00. Other work examines the stress field around a single simply connected stiff elastic 
inclusion subjected to a remote constant stress at infinity |23j and provides optimal lower bounds for 
the supremum of the maximum principal stress. The work presented in [5] provides an optimal lower 
bound on the supremum of the maximum principal stress for two-dimensional periodic composites 
consisting of a single simply connected elastically stiff inclusion inside the period cell. The recent 
work of [9] builds on the earlier work of [21 [15] and develops new lower bounds on the LP norm of 
the local stress and strain fields inside statistically isotropic two-phase elastic composites. However 
to date those bounds have been shown to be optimal for p = 2 see, [5] . Their optimality for p > 2 
remains to be seen. Optimal upper and lower bounds on the L 2 norm of local gradient fields are 
established using integral representation formulas in [16) . 

We conclude by providing the notation and summation conventions used in this article. Con- 
tractions of stress or strain fields a and e arc defined by a : e = Oijtij and \a\ 2 = a : a, where 
repeated indices indicate summation. Products of fourth order tensors C and strain tensors e are 
written as Ce and are given by [Ce]y = Cijutu] and products of stresses a with vectors v are given 
by [crvjj = &ijVj. The fourth order identity map on the space of stresses or strains is denoted by I 
and lijki = \/2(8ii~8ji + SuSjk)- The projection onto the hydrostatic part of <r(x) is denoted by ¥ H 
and is given explicitly by 

^H = \^kl and P*V(x) = ^)j. (1.1) 

The projection onto the deviatoric part of <j(x) is denoted by ¥ D and I = ¥ H + P D with P^P^ = 
f>HpD _ q rQjg t ensor product between two vectors u and v is the matrix u® v with elements[u <g) 
v]ij — UiVj. Last we denote the basis for the space of constant 3x3 symmetric strain tensors by e kl 
where e^„ = S mk S„i. 



2 Stress and strain fields inside heterogeneous thermoelastic 
media with imposed macroscopic loading 

Several distinct physical processes can generate prestress within heterogeneous media. In many 
cases it is generated by a mismatch between the coefficients of thermal expansion of the component 
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materials. To fix ideas we present the physical model associated with this situation. The tensors 
of thermal expansion inside each phase are given by Ai = h\I and A2 = h%I where I is the 3x3 
identity. The elastic properties of each component are specified by the elasticity tensors C\ and 
C2 respectively. In this treatment we consider heterogeneous elastically isotropic materials and the 
elasticity tensors of materials one and two are specified by 

d = 3ki¥ H + 2^P D , a = 1,2. (2.1) 

Without loss of generality we adopt the convention 

Ml > A*2- (2.2) 

The elastic displacement inside the composite is denoted by u and the associated strain tensor 
is denoted by e(u). The position dependent elastic tensor and thermal expansion tensor for the 
heterogeneous medium are denoted by C(x) and A(x) respectively where x denotes a point inside 
the medium. The domain containing the composite is given by a cube Q of unit side length. Here 
is supposed that Q is the period cell for an infinite clastic medium. In what follows the integral of 
a quantity q over the unit cube Q is denoted by (q) . 

A constant macroscopic stress a and uniform change in temperature AT is imposed upon the 
heterogeneous material. The local stress inside the heterogeneous medium is expressed as the sum 
of a periodic mean zero fluctuation a and ct, i.e., cr(x) = a + a(x) with (a) — 0. Elastic equilibrium 
inside each phase is given by: 

diva = 0. (2.3) 

The local elastic strain e(u) is related to the local stress through the constitutive law 

< r(x)=C(x)( e (u(x))-A(x)AT), (2.4) 

and the local elastic field is written in the form 

e(u) = e + e(u per ) (2.5) 

where u per is Q periodic taking the same values on opposite sides of the period cell and (e(u per )} = 0. 
Perfect contact between the component materials is assumed, thus both the displacement u and 
traction an are continuous across the two phase interface, i.e., 

u u = u, 2 , (2.6) 
a w n = a\ 2 n. (2.7) 

Here the subscripts indicate the side of the interface that the displacement and traction fields are 
evaluated on and n denotes the normal vector to the interface pointing from material one into 
material two. 

The effective "macroscopic" constitutive law for the heterogeneous medium is given by the con- 
stant effective elasticity tensor C e and effective thermal stress tensor H e that provide the linear 
relation between the imposed macroscopic stress a, uniform change in temperature AT, and the 
average strain e given by [18] , 

a = C e e + H e AT. (2.8) 

Here the components of C e are given by 

C&« = (Ci jm n(x)(e(<p kl ) mn + eti)), (2.9) 

where the fields are the periodic solutions of 

div{C{x){e{ip kl ) + e kl )) = inside each phase , (2.10) 

with the appropriate traction and continuity conditions along the two phase interface given by 

A = /V (2.11) 
Ci(e(^) + e fe '),n = C 2 (e(p H ) + e kl ) n. (2.12) 
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The effective thermal stress tensor H e is given by 

= <C(*)( e (^)-A(x))). (2.13) 

Where ip p is the periodic solution of 

div(C(x)(e(ip p ) - A(x))) = inside each phase , (2-14) 

with the traction and continuity conditions along the two phase interface given by 

A = V P | 2 , (2.15) 
d (e(^) - Ai) u n = C 2 (e(^)-A 2 ) |2 n. (2.16) 

From linearity it follows that the local fluctuating strain field is given by the sum 

e(u per ) = t{ip kl )l kl + e(/)AT. (2.17) 

We write tp e — ip kl e kl and from linearity one has the expression for C e e given by 

C e e=(C(x)( e (^) + e->. (2.18) 

In this article the imposed mechanical stress is given by a constant hydrostatic stress 

a = <7 I, (2.19) 

where <7o can assume any value in — oo < <7q < oo. We introduce a method for obtaining optimal 
lower bounds on the higher moments of the hydrostatic component of the local stress inside the 
composite when it is subjected to an imposed hydrostatic load ctqI and temperature change AT. 
Here no restriction is placed on AT. The volume fractions of materials one and two are denoted by 
9\ and 9 2 and the average of a quantity q over material one is denoted by (g)i and over material 
two by {q)2- In the following section we present optimal lower bounds on the following moments of 
the local hydrostatic stress P H a(x) over the domain occupied by each material given by 

(|P ff fr(x)Hi /p and (|P ff a(x)| p )2 /p , (2.20) 
for 1 < p < oo, as well as for the maximum local hydrostatic stress over the whole composite domain 

max{|P H CT(x)|} . (2.21) 

It is pointed out that the case corresponding to p = oo in (|2.20l) corresponds to lower bounds on 
the maximum local stress over each phase 

max ||P H cr(x)|}, max {|P ff cr(x)|}, (2.22) 

X in material 1 L J X in material 2 L J 

The lower bounds are given in terms of the volume fractions 6\ and 02, as well as the bulk and 
shear moduli k±, k2, Mi, fJ-2, and the coefficients of thermal expansion hi and \i2- The lower bounds 
are described by the following characteristic combinations of these parameters given by: 

= fci(fc 2 + |M2) 

1 k l k2 + {k 1 6 1 + k 2 6 2 )^2 ' 



L = 3 (2 241 

fclfc 2 + (Ml+fc202)f/V 

ki{k 2 + f/^i) 
hk 2 + (ki0i + k 2 e 2 )^i 



k\ (fc 2 + qWi) 



fc 2 (fci + |/i 2 ) 

M 2 = ^ ' . , 2.26 

k 1 k2 + (k 1 9 1 +k 2 92)l^2 ^ ' 
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D = AT( v -!-), (2.27) 

fc 2 - fci 

and 

^ = ^(1-T7TM7)- ( 2 - 28 ) 

2 

For elastically well-ordered materials, ki > k2 one has L\ > 1 > L 2 , and Mi > 1 > M2; for the 
non well-ordered case > fci, one has L 2 > 1 > Li and M 2 > 1 > Mi. 

The lower bounds are shown to be obtained by the local fields inside the coated sphere assemblages 
introduced in [7]. The lower bounds presented here include the effects of thermal stresses due to 
thermal loads and reduce to the optimal bounds reported in [15) when AT = 0. 

3 Lower bounds on local stress for elastically well-ordered 
t her mo elastic composite media 

In this section it is assumed that the materials inside the heterogeneous medium are elastically 
well-ordered, i.e., /ii > fi2 and ki > K2- We present lower bounds that are optimal for the full 
range of imposed hydrostatic stresses, i.e., —00 < cto < 00 as well as for unrestricted choices of AT. 
The configurations that attain the bounds are given by the coated sphere assemblages [7]. To fix 
ideas we describe the coated sphere assemblage made from a core of material one with a coating of 
material two. We first fill the cube Q with an assemblage of spheres with sizes ranging down to the 
infinitesimal. Inside each sphere one places a smaller concentric sphere filled with "core" material 
one and the surrounding coating is filled with material two. The volume fractions of material one 
and two are taken to be the same for all of the coated spheres. 

In what follows we list the lower bounds for the well ordered case. These bounds are derived 
in Section 5. Their optimality follows from explicit formulas for the moments of the local fields 
inside the coated sphere assemblage, these are discussed and presented in Section 5. The first set 
of bounds apply to all moments (|P ff cr| p ).|/ p for 1 < p < 00. We suppose that h-2 > hi fix AT and 
list the bounds as a function of the imposed macroscopic stress a®. The bounds are displayed in 
the following table where the optimal microstructures are given by the coated spheres construction. 
The coating and core phase of the optimal configuration is listed in the table below. 



Range 


LowerBound 


Optimal 
microstructure 


—00 < do < D 


{\¥ H a\P )l /p > V3[(D - a Q )L2 - D] 


Core material 2 and 
coating material 1 


D<a <D(l-^- 2 ) 


<|P"a|P) 2 /p > V3[(D - a )M 2 - D] 


Core material 1 and 
coating material 2 


D(l-^)<a <D(l-±) 


{\f H a\P)\ lP > 


Optimality undeter- 
mined 


D(l - < (To < 00 


(\F H a\P)l /p > s/3[(a - D)L 2 + D] 


Core material 2 and 
coating material 1 



Next we suppose that hi > ft-2 and present optimal lower bounds on (\P H a\ p )2 , for 1 < p < 00. 
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The bounds and associated optimal microstructures are given in the following table. 



Range 


LowerBound 


Optimal 
microstructure 


-co < a < D(l - 


(\¥ H a\P )\ /p > VS[(D - a Q )L 2 - D] 


Core material 2 and 
coating material 1 


D(l-±)<* <D(l-±) 


(\¥ H a\P) 1 2 /p > 


Optimality 
undetermined 


D(l -^)<a <D 


(\¥ H a\P )l /p > V3[(a - D)M 2 + D] 


Core material 1 and 
coating material 2 


D < ctq < 00 


(\¥ H a\P )l /p > V3[(a - D)L 2 + D] 


Core material 2 and 
coating material 1 


Lower bounds and the associated optimal microstructures for all moments (|P' Ef a-| p } 1 ^ p for 1 
< 00 for the case h 2 > hi is given in the following table. 


Range 


LowerBound 


Optimal 
microstructure 


—00 < <to < D 


(\¥ H a\P)l /p > V3[(D - ao)^ - D] 


Core material 1 and 
coating material 2 


D < a < D(l - -fa) 


(\¥ H a\P)\ /p > VS[(D-a )M 1 -D] 


Core material 2 and 
coating material 1 


D(l-^)<a <D(l-±) 


(\¥ H a\P)l /p > 


Optimality 
undetermined 


D(l - ±) < (T < OO 


(\¥ H a\P)\ /p > V3[(a -D)L 1 + D} 


Core material 1 and 
coating material 2 



Lower bounds and the associated optimal microstructures for all moments (\F H cr| p ) 1 for 1 < p < 00 
for the case hi > h 2 is given in the following table. 



Range 


LowerBound 


Optimal 
microstructure 


-00 <<7 < D(l - -i-) 


(\¥ H a\P)l /p > VS[(D - a )Li - D] 


Core material 1 and 
coating material 2 


D(l-±)<v <D{l-±) 


(\¥ H a\P) 1 i /p > 


Optimality 
undetermined 


D(l -^)<<J <D 


(\¥ H a\P)\ /p > V3[(a -D)Mi+D] 


Core material 2 and 
coating material 1 


D < ctq < 00 


(\¥ H a\P)l /p > VS[(a -D)Li+D] 


Core material 1 and 
coating material 2 



Next we display lower bounds on max xin Q {|P ff (r(x)|}. We start with the case h 2 > hi and 
the lower bounds and optimal geometries are given in the following table. The phase in which the 
maximum is attained is denoted with an asterisk. 



Range 


LowerBound 


Optimal 
microstructure 


—00 < (To < D 


max XlnQ {|P W (r(x)|} > V3[(D - a )Li - D] 


Core material 1* and 
coating material 2 


D<<r <F 


max XlnQ {|P ff (r(x)|} > V3[(D - a )M 2 - D] 


Core material 1 and 
coating material 2* 


F < (T < 00 


max XlnQ {|P W (r(x)|} > V3[(a - D)Li + D] 


Core material 1* and 
coating material 2 
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Lower bounds for and optimal microgepmetries for the case hi > h 2 are given in the following table. 
The phase in which the maximum is attained is denoted with an asterisk. 



Range 


LowerBound 


Optimal 
microstructure 


— OO < <7q < F 


max xin q {|P ff (r(x)|} > \/3[(D — uq)Li — D] 


Core material 1* and 
coating material 2 


F <a <D 


max xm q {|P H a(x)|} > V3[(o- - D)M 2 + D] 


Core material 1 and 
coating material 2* 


D < ctq < 00 


max xinQ {|P^(7(x)|} > V3[(<i - D)L 1 + D] 


Core material 1* and 
coating material 2 



4 Lower bounds on local stress for non well-ordered ther- 
moelastic composite media 

In this section it is assumed that the materials inside the heterogeneous medium are elastically non 
well-ordered, i.e., \i\ > fi 2 and k 2 > n±. We fix AT and present lower bounds that are optimal for 
the full range of imposed hydrostatic stresses, i.e., —00 < cto < °°- The configurations that attain 
the bounds for the non well-ordered case are also given by the coated sphere assemblages [7]. 

In what follows we list the lower bounds for the non well-ordered case. These bounds are derived 
in Section 5. Their optimality follows from explicit formulas for the moments of the local fields 
inside the coated sphere assemblage, these are discussed and presented in Section 5. The first set 

of bounds apply to all moments (\F H a\ p )2^ P for 1 < p < 00. We suppose that h 2 > hi and list 
the bounds as a function of the imposed macroscopic stress cto. The bounds are displayed in the 
following table where the optimal microstructures are given by the coated spheres construction. The 
coating and core phase of the optimal coated sphere configuration is listed in the table below. 



Range 


LowerBound 


Optimal 
microstructure 


-00 < a < D(l - jfe) 


{\¥ H a\P)l /p > V3[(D - a )M 2 - D] 


Core material 1 and 
coating material 2 


D(l-fa)<a Q <D{l-£) 


(\¥ H a\P}l /p > 


Optimality 
undetermined 


D(l ~^-)<<j <D 


{\F H a\P )\ /p > V3[(a ~ D)L 2 + D] 


Core material 2 and 
coating material 1 


D < (To < 00 


{\F H a\P )\ lp > V3{(a - D)M 2 + D] 


Core material 1 and 
coating material 2 



Next we suppose that hi > h 2 and present optimal lower bounds on (\¥ H a\ p )2 , for 1 < p < 00. 
The bounds and associated optimal microstructures are given in the following table. 



Range 


LowerBound 


Optimal 
microstructure 


—00 < cto < D 


{\F H a\P )\ /p > V3[(D - ct )M 2 - D] 


Core material 1 and 
coating material 2 


D < a a < D(l - ±) 


{\F H a\P )\ lp > V3[(£> - ct )L 2 - D] 


Core material 2 and 
coating material 1 


D{l-±)<o- <D(l-^) 


(\F H a\ p ) 1 2 /p > 


Optimality 
undetermined 


D(l -^)<a <oo 


(|P"ct|^ /p > s/3[(a - D)M 2 + D] 


Core material 1 and 
coating material 2 
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Lower bounds and the associated optimal microstructures for all moments (\F H a\ p } 1 for 1 < 
p < oo for the case h 2 > hi is given in the following table. 



Range 


LowerBound 


Optimal 
microstructure 


-OO < <T < D(l - -fc) 


(\¥ H a\P)l /p > V3[(D-a )M 1 -D] 


Core material 2 and 
coating material 1 


D(l - < tro < D(l - 


(\¥ H a\P)l /p > 


Optimality 
undetermined 


D(l - A-) < <t < D 


(|P H a|f)l /p > x/3[(<t -£>)Li+£>] 


Core material 1 and 
coating material 2 


D < <Tq < 00 


(|pff CT |p)VP > V3[(a -D)M 1 +D] 


Core material 2 and 
coating material 1 



Lower bounds and the associated optimal microstructures for all moments (\F H a\ p )\^ p for 1 < p < 00 
for the case hi > hi is given in the following table. 



Range 


LowerBound 


Optimal 
microstructure 


—00 < a < D 


(\¥ H a\P)l /p > VS[(D - a Q )Mi - D] 


Core material 2 and 
coating material 1 


D < <j < D(l - 


{\¥ H a\P)\ /p > VS[(D- a Q )Li - D] 


Core material 1 and 
coating material 2 


D(l-±)<a <D(l-±) 


{\W H a\v)\ lp > 


Optimality 
undetermined 




(\¥ H a\P)\ /p > VS[(a - D)Mi + D] 


Core material 2 and 
coating material 1 



Next we display lower bounds on maxxm q 

{|P ff a(x)|}. We start with the case hi > h\ and 
the lower bounds and optimal geometries are given in the following table. The phase in which the 
maximum is attained is denoted by an asterisk. 



Range 


LowerBound 


Optimal 
microstructure 


— OO < <T{) < F 


max xin q {\¥ H a(x)\} > V3[(D - a )M 2 - D] 


Core material 1 and 
coating material 2* 


F <<J <D 


max Xln q {|P"a(x)|} > VS[(a - D)L X + D] 


Core material 1* and 
coating material 2 


D < <To < 00 


max xin q {|P w ct(x)|} > V3[(a - D)M 2 + D] 


Core material 1 and 
coating material 2* 



Lower bounds for and optimal microgcpmetries for the case hi > hi are given in the following table. 
The phase in which the maximum is attained is denoted by an asterisk. 



Range 


LowerBound 


Optimal 
microstructure 


—00 < (To < D 


max XlnQ {|P"(t(x)|} > y/3[(D - a )M 2 - D] 


Core material 1 and 
coating material 2* 


D<cr <F 


max XlnQ {|P ff (r(x)|} > VS[{D - a )Li - D] 


Core material 1* and 
coating material 2 


F < (To < 00 


max Xin q {|P H (t(x)|} > V3[{a - D)M 2 + D] 


Core material 1 and 
coating material 2* 
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5 Derivation of the lower bounds on < |P#cr| p >f 



In this section we outline the methodology for proving optimal lower bounds. The bounds are 
derived using duality relations. We use the following duality relation posed over the space of square 
integrable symmetric matrix fields 77 that holds for p > 1 given by 



h\F H a\t>) i =sup\{F H a:r ] } i -± J {\r l f) i \, fori = 1,2, (5.1) 
v v I P ) 



where p' is the conjugate exponent to p given by p' — ^zy- This relation follows immediately from 
standard duality relations see, Restricting 77 to the set of all constant matrices and taking the 
supremum delivers the basic bounds: 

(\¥ H a\ p ) l > |(P H CT)i| p , p> 1 and for i= 1,2. (5.2) 

We point out that equality holds in (|5.2[) if and only if P H a is identically constant inside the i th 
material. In what follows we outline the method for obtaining bounds on the moments in material two 
noting that the identical procedure delivers bounds on the moments in material one. We introduce 
the indicator function of material one xi taking the value 1 in material one and zero outside. The 
indicator function corresponding to material two is denoted by X2 and X2 = 1 — Xi- To proceed we 
rewrite the right hand side of (15. 2| in terms of the effective elastic properties and thermal expansion 
coefficient. To do this we use the following identity given by 

q 7^ 

tr < X2cr >= 2 —(a - AWCT 1 / : / + k 1 AT(C e r 1 H e : / + fciAT < A >: /). (5.3) 

fc 2 - ki 

This identity is obtained in the following way. Taking averages on both sides of (|2.4p gives 
< a >= (C(x)(e(u(x)) - AT\(x))). On writing C(x) = C x + (C 2 - Ci)x 2 (x) we see that 

< a >= C x (e- < A >) + (C 2 - Ci)C^ < X20 > ■ (5.4) 

From (|2~8|) we see that e = (C e ) _1 (< a > -ATH e ) and for (a) = a I we obtain 

< X20 >= C 2 (C 2 - CO-Vo/ - Ci((C e )-V I - AT(C e )~ 1 H e — AT < A >)). (5.5) 



The identity (|5.3p now follows by applying the hydrostatic projection P# to both sides of (|5.5j) . 
We now derive the lower bound. Applying the basic bound (|5.2|) to (|P//cr| p ) 2 and ()5.3|) gives 



tr(( X 2<r))Y 



>l ^ >21 ^' V3, 2 

= 3 p/2 g" p | , k \ \ p \kia ( i - (C 6 )- 1 / : i) + kiAT(C e )~ 1 H e I : / + fciAT(A) : 7| p .(5.6) 
k 2 - ki \ ki J 

We note that equality holds in (|5.6p when Pj^cr is constant inside material two. 

We now employ an exact relation that relates the contraction (C e )~ 1 H e : / involving the effective 
thermal stress tensor H e to the quantity (C e ) -1 / : /. The exact relation used here is given by 

3(/i 2 - hiXC 6 )- 1 / : I + 3(^ - 4*) 
(C 6 )" 1 ^ 6 :/ = — 1!—-^ j 1^2 (5.7) 

fel A'2 

This exact relation is a direct consequence of the exact relation developed by [22] for the effective 
thermal expansion tensor a e = —(C e )~ 1 H e . 

Substitution of (|5.7|) into (|5.6p and algebraic manipulation gives 

(|P^| P ) 2 /P > V3|((ro- r>)X + £>|, (5.8) 
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where 

*-*=^(ir< crlJ:J )- (5 ' 9) 

As before we point out that equality holds in (|5.8[) when Pjycr is identically constant inside material 
two. 

Identical arguments deliver the lower bound on the moments over material two given by 



p\i/p 



i > VS\(<t -D)Y + D\, (5.10) 



where 



' i 



where (|5.10|) holds with equality when F H <r is a constant in material one. 

The variables X and Y are constrained to lie within intervals set by bounds on the contraction 
(C e ) -1 / : /. These bounds follow immediately from the work of Kantor and Bergman [TT] and are 
given by 

(K+ s T 1 <(C e r 1 I--I<(K HS )-\ (5.12) 



,,,,, ,, r,™ + U„ t 



where K HS and are the Hashin and Shtrikman bulk modulus bounds [5] given by 



k+ s = Ml + - ( //f;;" i} 4 ) (5.i3) 

fcl0 2 + fc2#l + 3MI 

and 

K HS = k 1 e 1 + k 2 6 2 - {————— — )• 5.14) 

ki6 2 + k 2 0i + ffi 2 

These bounds hold both for elastically well-ordered materials and elastically non well-ordered ma- 
terials. When the materials are well-ordered (|5.12p implies that X and Y lie in the intervals 

L 2 < X < M 2 , (5.15) 

L x < Y <M X . (5.16) 

while for non well-ordered materials 

M 2 < X < L 2 , (5.17) 

Mi < Y <L X . (5.18) 

A straightforward calculation in Section [6] shows that the hydrostatic component of the local 
stress is constant inside each phase of the coated sphere construction. Hence ()5.8|) and ()5.10|) hold 
with equality for the coated spheres construction and we obtain explicit formulas for the moments 
of the hydrostatic component of the local stresses for these composites. For coated spheres with 
core phase 1 and coating phase 2, (C 6 )" 1 / : I = (K^g) -1 and substitution of (|5.14[) into (|5.9|) and 
(|5.1ip together with (|5.8[) and (|5.10p shows that the moments are given by 

{\¥ H a\ p ) 2 /p = Vs\(a a -D)M 2 +D\, and (5.19) 

(|P^| P )1 /P = V3 |(tro - D)Li + D\ . (5.20) 

For coated spheres with core phase 2 and coating phase 1, (C e ) -1 / : I = (K^g)" 1 and substitution 
of (|5.13p into (|5.9p and (|5.1ip together with (|5.8p and (|5.10p shows that the moments are given by 



{\F H cx\ p ) 1 2 /p = V3\(a -D)L 2 + D\, and (5.21 

(|Pk0| p )i /p = V3|(tr - D Wi + D \ ■ (5-22 
We collect results and state the lower bounds and indicate when they are optimal. 
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Theorem 5.1. Bounds for well-ordered composites, k\ > k 2 . For 1 < p < oo, any choice of AT, 
and — oo < (To < oo the lower bounds are given by the following formulas. 

(|PW>2 /P > T mm {V3\(a -D)X + D\\, (5.23) 

L 2 <X<M 2 ( J 

and when the minimum is realized for X = L 2 the bound is attained by the fields inside the core 
phase of a coated sphere construction with core material 2 and coating 1; when the minimum is 
realized for X = M 2 the bound is attained by the fields inside the coating phase of a coated sphere 
construction with core material 1 and coating 2. 



<|P^nl /p > Li< mjn Mi {V3 |(a - D)Y + D\) , (5.24) 

and when the minimum is realized for Y = L\ the bound is attained by the fields inside the core 
phase of a coated sphere construction with core material 1 and coating 2 and when the minimum is 
realized for Y = M\ the bound is attained by the fields inside the coating phase of a coated sphere 
construction with core material 2 and coating 1. 

Theorem 5.2. Bounds for non well-ordered composites, k 2 > k\. For 1 < p < oo, any choice of 
AT, and — oo < cr < oo the lower bounds are given by the following formulas. 

(\PHO-\nl /p > M2 min <L2 {V3 \(a - D)X + D\} , (5.25) 

and when the minimum is realized for X = L 2 the bound is attained by the fields inside the core 
phase of a coated sphere construction with core material 2 and coating 1; when the minimum is 
realized for X = M 2 the bound is attained by the fields inside the coating phase of a coated sphere 
construction with core material 1 and coating 2. 



(\^ho-\ p )\ /p > min {V3|(a - D)Y + D\\ , (5.26) 

and when the minimum is realized for Y = L\ the bound is attained by the fields inside the core 
phase of a coated sphere construction with core material 1 and coating 2 and when the minimum is 
realized for Y = M\ the bound is attained by the fields inside the coating phase of a coated sphere 
construction with core material 2 and coating 1. 

These bounds are stated explicitly in the first four tables of Sections 3 and 4. 

We conclude by outlining the steps behind the derivation of the lower bounds on the maximum 
values of the local fields inside thermally stressed composites. For the well-ordered case we use the 
simple lower bound given by 

max {|P H cr(x)} > ui&k{A,B}. (5.27) 

X in Q 

where 

A = L2 nAn M2 {V3\(vo-D)X + D\}, 

B = min \V3\(a -D)Y + D\\. (5.28) 

Li<F<Mi I J 

For the non well-ordered case we use 

max {|P H o-(x)} > max{C,L>}. (5.29) 

X in Q 

where 

c = mM\J v ~ 3 ^°- d)x + d] }> 

D = min \VS\(a Q -D)Y + D\). (5.30) 

M 1 <Y<L 1 L J 
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The bounds given in the last two tables presented in Sections 3 and 4 follow from straight forward 
but tedious calculation of the explicit formulas corresponding to (|5.27[) and (|5.29[) . A delicate but 
straight forward computation shows that these lower bounds are attained by the fields inside the 
coated sphere assemblage. 



6 Local stress and strain fields inside thermally stressed coated 
sphere geometries 

In this section we summarize the properties of local fields inside the coated sphere assemblage in 
the presence of thermal stress due to a mismatch in the coefficients of thermal expansion. From 
linearity the local stress can be split into the sum of two components; one component arising from 
imposed mechanical stress and a second component associated with thermal stress. It is known that 
the local stress due to an imposed hydrostatic stress has constant hydrostatic part inside each phase, 
this follows from explicit solution see for example [18] . Here we display the explicit solution for the 
local stress due to mismatch in the coefficients of thermal expansion and show that it has a constant 
hydrostatic component inside each phase. From this we conclude that the total local stress inside 
the coated sphere assemblage has a constant hydrostatic component inside each phase. 

We solve for the stress inside a prototypical coated sphere composed of a spherical core of material 
two with radius a, surrounded by a concentric shell of material one with an outer radius b. The ratio 
(a/b) 3 is fixed and equal to the inclusion volume fraction 8 2 . Here the the coefficients of thermal 
expansion for the core and coating are given by hi and h 2 respectively. The local elastic displacement 
(p satisfies the equations of elastic equilibrium are given by: 



div(C 2 (e(£) - h 2 I)) = 
cfo(CiO#) - hil)) = 
Ci(e(0)- W)n|i = C7 2 (e(0) 
(p is continuous 

(f> = 



h 2 I)n\, 



< r < a, 
a < r < b, 

continuity of traction at r = a, 

on < r < b, 

on the boundary r = b. 



We assume a general form of the solution given by 



Cr < r < a, 

<p = { Ar + Bf; a<r < b, 
r > b, 

where r = |r|, n = r/r and A, B, C are unknowns. The corresponding strain field e((p) is given by 



CSij < r < a, 

(e(<p))ij = ^(<Pi,j + <Pj,i) = I A5 l3 + ^(Sij - in^j) a <r <b, 

r > b. 



(6.1) 



On applying the continuity of displacement and the traction at the interface we find that 

36> 2 (fci/ii - k 2 h 2 ) 



A 



B 



C 



3M 2 +4/xi + 3fc 2 (l-0 2 )' 

— 3a 3 (k\hi — k 2 h 2 ) 
3k 1 6 2 + 4^H+3k 2 (l-6 2 ) , 

-3(1 - e^jhhi - k 2 h 2 ) 
3fci6»2+4^i+3fc 2 (l-6»2)' 
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Computation of the radical component of the stress at r = b gives 



Ci(e(<p) - h x I)n = H*n 



(6.2) 



where H* 



H* = 



36> 2 (3fci+4^i)(fcifei-fc 2 fe2) 



I — 'Skihil. 



(6.3) 



In this way we have constructed a solution (p for the elastic field inside every coated sphere in 
the assemblage. We now define (p p on the whole domain Q to be given by <p inside each coated 
sphere and zero outside. It easily follows on integrating by parts using ()6.2j) together with the 
fact that if p vanishes on the boundary of each coated sphere that ip p is the weak solution [4] of 
div(C(e(ip p — A)) = over the full domain Q, i.e., 



for every periodic test function <j). Equation (|6.1[) implies that hydrostatic component of stress is 
constant inside each phase. 

Last we show that the effective thermal stress H e for the coated sphere assemblage is given by 
H*. Inside each coated sphere Si, i = 1,2, ... we consider a = C(e(ip p ) — A) and integrate by parts 
and apply (16. 2| to find that 



where ds is an element of surface area on the outer surface of the coated sphere dSi and \Si\ is the 
volume of Si. Substitution of (|6.5j) into (|2.13j) gives the required identity 
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